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We study the electromagnetic form factors, decay constants and charge radii of the pion and 
kaon within the framework of light-front field theory formalism where we use an ansatz for the 
quark-meson interaction bound-state function which is symmetric under exchange of quark and 
antiquark momentum. The above mentioned observables are evaluated for the -I- component of the 
electromagnetic current, J'^, in the Breit frame. We also check the invariance of these observables 
in other frames, whereby both the valance and the non-valence contributions have to be taken into 
account, and study the sensitivity of the electromagnetic form factors and charge radius to the 
model’s parameters; namely, the quark masses, ms, and the regulator mass, mu. It is 

found that after a fine tuning of the regulator mass, i.e. mu = 0.6 GeV, the model is suitable to ht 
the available experimental data within the theoretical uncertainties of both the pion and kaon. 

PACS numbers: 


I. INTRODUCTION 

The theory of strong interactions. Quantum Chromo¬ 
dynamics (QCD), has been the object of theoretical and 
experimental scrutiny for four decades now and is per- 
turbatively well defined. Indeed, at large momentum 
transfer, perturbative calculations successfully describe a 
wealth of subatomic phenomena. However, QCD is also 
a theory whose elementary excitations are confined and 
even in the far infrared below the typical hadronic scale, 
E Ki 1 GeV, QCD does not seem to break down [1-3]. 
Thus, nonperturbatively QCD may well be rigorously de¬ 
fined but a full solution is yet out of reach: there exists no 
simple Schrodinger picture of a many-body Hamiltonian 
as in quantum mechanics [4]. This is due to the intrin¬ 
sic nonperturbative nature of quark-antiquark pair cre¬ 
ation and annihilation in a relativistic quantum field the¬ 
ory, which entails non-conservation of particle number. 
In hadron physics, this situation leads to seek computa¬ 
tional approaches beyond perturbative QCD and many 
models or effective theories have been proposed to tackle 
QCD in the nonperturbative regime with the aim to de¬ 
scribe hadron phenomenology. 

In this context, one possibility to develop a nonpertur¬ 
bative covariant framework is the light-front field theory 
formalism proposed by Dirac in 1949 [5]. In this ap¬ 
proach, the hadronic bound states are described by wave 
functions in the light-front space-time hyper surface, de- 
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fined by the coordinate x'^ = -\- = 0, and due to 

the stability of the Fock-state decomposition these wave 
functions are covariant under kinematic boosts [4, 6]. 

On the other hand, in understanding the dynami¬ 
cal properties of nonperturbative QCD, the light pseu¬ 
doscalar mesons and in particular the pion play a crucial 
role. Remarkably, the latter is a bound state of mas¬ 
sive antiquark-quark pairs as well as the almost massless 
Goldstone boson associated with chiral symmetry break¬ 
ing. There have been many studies of their static proper¬ 
ties [7-12] and their dynamical properties have also been 
investigated theoretically [13-17, 24-62] and experimen¬ 
tally [64-74]. 

Taking advantage of the simple structure of the Fock 
space and the vacuum in light-front quantization, vari¬ 
ous hadronic properties of bound states, such as decay 
constants and electromagnetic form factors of the pion, 
kaon and nucleon, have been calculated [16, 17, 52-62] 
and successfully compared with their experimental val¬ 
ues [64-74]. Since the light front component, J+, has 
been successfully employed to calculate electromagnetic 
form factors [7, 31, 63, 75-79], the light-front approach 
also offers a theoretical framework to extract from them 
useful information on the valence and non-valence com¬ 
ponents of the meson’s wave function. 

In the present simultaneous study of electromagnetic 
form factors, charge radii and decay constants, we adopt 
the light-front field theory formalism of Refs. [7, 15] 
wherein the Bethe-Salpeter amplitude of the qq bound 
states was modeled for two different momentum con¬ 
stellations, namely a symmetric [17] and nonsymmet- 
ric vertex model [16]. Here, the vertex refers to the 
qq pair coupling to the pseudoscalar meson in an effec¬ 
tive Lagrangian. Using a nonsymmetric vertex model, 
E. O. Silva et. al. [58] calculated the aforementioned 
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pion and kaon observables which are in agreement with 
experimental data. However, a momentum distribution 
of the meson that is symmetric under the exchange of the 
quark and antiquark momenta is more realistic and such 
a model for the Bethe-Salpeter amplitude should improve 
the description of or at least equally well reproduce all 
observables presented in Ref. [58]. Thus, we here use the 
same component, J"*", of the light-front electromagnetic 
current, though with a symmetric momentum descrip¬ 
tion of quark-meson bound-state vertex. It is important 
to recall here that the choice of J’*' with the Drell-Yan 
condition g'*" = 0 guarantees that pair-term contribu¬ 
tions (non-valence terms) vanish [16, 17]. On the other 
hand, to preserve rotational symmetry, the pair contri¬ 
bution must be included [18-23]. Consequently, we here 
employ both the valence and non-valence contributions 
considering the case q'^ ^ 0. 

This paper is organized as follows: Sec. II serves to 
summarize the general framework where subsequently 
the different physical observables, namely the electro¬ 
magnetic form factors, decay constant and charge radii 
are discussed in turn. In Sec. Ill, we present our numer¬ 
ical results and analyze the observables’ dependence on 
variation of the model parameters. In the last section, 
we give our conclusions. 


II. THE MODEL 

In this section, we briefly summarize the model and the 
computational tools of the light-front formalism required 
to investigate the pseudoscalar meson’s electromagnetic 
form factors, charge radii and the decay constants. Our 
approach is based on similar earlier work [16, 17], where 
the following effective Lagrangian for the qq bound state 
was employed: 

Ces = -igcj)-q-f^rq , (1) 

where g = too-//o- is the coupling constant, ttiq- and 
/o- denote the mass and decay constant of a pseudoscalar 
meson, respectively, and cj) represents the scalar field. We 
make a symmetric ansatz for the gg-meson vertex which 
describes the bound state, 

A(fc, P) = C [(fc^ — m|j -I- ie)~^ 

+ {{P-kf+ ( 2 ) 

where it is clear that the A(k, P) is symmetric under the 
exchange of the quark and antiquark momenta, k and 
{P — k)\ P is the total momentum of the meson and C 
a normalization constant. In the following, we discuss 
the analytic light-front formulation of the electromag¬ 
netic form factor, charge radius and decay constant. 


A. Electromagnetic form factors 


rent, = e^gy^g, is sandwiched between the initial and 
final bound states of the same meson: 


P^PoY“(g2) = (Mo-(p')I4|Mo-(p)), (3) 

where Mq- = 7r+, Ar+, P^ = {p + p')^ and g^ = ip-p'Y 
is the square of the momentum transfer. 

The electromagnetic form factor in the impulse approx¬ 
imation is obtained from triangle diagrams, each of which 
contains one spectator quark. In this approximation, 
the covariant electromagnetic current of a pseudoscalar 
meson, J^, that enters Eq. (3) can be written as fol¬ 
lows [80, 81]: 


= N 


d'^k 

1 


Tr 


1 


1 


-7 


7m 


’ - p-m„ 


— niq + ie ' ^ f — iriq + le 

— 7 ® A{k,p')A{k,p) -k [g o g], (4) 


with the normalization. 


-2ieqml„Nc „ m„ + mq 

N = - 72^ - : ™ 0 - : - 5 - , 

h- ^ 

where Nc = 3 is the number of colors and /q- the pseu¬ 
doscalar weak decay constant. 

The light front-form variables are. 


k^ = k^ + k^, k =k^ — k^, k± = {k^,k'^) 
g"*" = \/—g^ sin a, q^ = \/—q^ cos a, qy = 0 

g2 = q+q- _ (gl)2 . (5) 


Here, we use the Drell-Yan condition, g+ = 0, in the 
Breit frame, which implies a = 0. However, our results 
are frame invariant, i.e. invariant for a ^ 0 where both 
the valence and the non-valence contributions become 
important. After introducing the front-form variables in 
Eq. (5) and using 7 + = 7 ° -|- 7 ^ to obtain the compo¬ 
nent of the current in Eq. (4), the electromagnetic form 
factor becomes: 

I ^‘‘k+dk-Tvp+] 

X r(A:+,p+,p'+) -k [g O g] , (6) 


In Eq. ( 6 ), the trace in light-front coordinates is. 


iTr[e>+] = ifc+g^-k (fc+-p+-p'+)(A:i -/c+fc") 

— — {p''^k± ■ p± + p~^k± ■ p'j_) 

— fc+(2fc^ -km? — 2mqmq) — (p+ + p''^)mqmq, 

( 7 ) 


and 

r{k+,p+,p’+) = 


A(fe+,p+)A(fc+,p'+) 


— rriq + ze)((p — fc)^ — m^ -k ie) 
1 


( 8 ) 


The covariant electromagnetic form factor, P®P is de¬ 
fined by a matrix element where the electromagnetic cur- 


X 


{{p' — fc )2 — mq+ ie) 
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where = k~^{k —k^^) where k^^ is the on-energy-shell Differentiation with respect to of the above equation 
value of the corresponding momentum given by, yields the charge radius, 


kon = 


kl 


k+ 


(9) 


In terms of light-front variables, the bound-state function 
in Eq. (2) becomes. 


A(fc+,p+) 


C 


{p+ - k+){p- - k- - 


p+—k+ 


) 


k+{k- 


c _ 

—it 

fc+ 


( 10 ) 


Collecting all ingredients from Eqs. (7)-(10), we insert 
them in Eq. (6) and after k~ energy integration (see ap¬ 
pendix of Ref. [17]) with x = ^ the electromagnetic 
form factor can be rewritten, 


AA r rl'^k , rlr 

Fr = J ^^^^<^*{x,k^Mx,k^)e{x)e{l-x) 

X - {p'^k± ■ PJ_ +p+k± ■ p'j_)] 

( 11 ) 

where AA = and 


d>(x, fe±,p+,p±) = 


1 


(1 - x)(ml_ - m|j)) 

1 


x(mg_ — A4^(m^, m|)) J mg_ — Ad^(m2, ^^2^ ’ 


+ [q^q] 


( 12 ) 


with 




2/_2 2 \ 

rrifc) = 




(p-k) 


■ml 


(1-a:) 


P\ 


Note that the appearance of the second term of Eq. (12) is 
due to the symmetric character of the meson-quarks ver¬ 
tex, which is absent in Refs. [16, 58] where the authors 
consider a nonsymmetric behavior of the vertex function. 
As it was shown [16, 17, 25], to preserve general covari¬ 
ance the non-valence contribution is mandatory. Thus, in 
Eq. (11) the step functions 9{x) and 0(1 — x) delimit the 
integration interval, 0 < fc < p"*", of the valence contribu¬ 
tion, whereas the interval, p'^ < k'^ < p''^, corresponds 
to the non-valence contributions to electromagnetic cur¬ 
rent [17]. 


B. Charge radius and decay constant 

The mean-square electric charge radius of a meson is a 
relevant quantity and correlated with the electromagnetic 
form factor, 

Fr{q^)c^l-^{rl)q^ . (13) 


jpem 

/j.2 \ _ “-^0- 

^ 0 -^" dq^ 


(14) 


9^=0 


A relevant observable and also our main constraint on 
the model’s parameters is given by the weak decay con¬ 
stant, /o-. The decay constant of a qq, bound state can 
be found from the following matrix element of the par¬ 
tially conserved axial-vector current: 


{0\A^^\0-)=^fo-p^^ . (15) 


where = ( 77 ^ 7 ^ 9 , is the axialvector current. The weak 
decay constant is given by. 


fo- 


fo- 


d‘^k 

(M 


rTr 


ie 


5 

(fc — p)2 — m? -I- ie 


A{k,p) 


(16) 


We make use of the -I- component of the axialvector cur¬ 
rent A'^ and after integration over k~, one obtains the 
decay constant in terms of the valence component of the 
model: 


fo- = y [4xTOq-|-4mg(l - x)]4>(x, fc_L,m^, 0 ). 

(17) 

III. NUMERICAL RESULTS AND DISCUSSION 

The model introduced in Sec. II contains three free pa¬ 
rameters, namely, the regulator mass, m^, and the two 
constituent quark masses, = md and mg, where the 
strange quark mass is taken from the study in Ref. [24]. 

The main focus of this study is to constrain the param¬ 
eters of a more realistic bound-state ansatz to accom¬ 
modate the available experimental data on the pion and 
kaon elastic form factors, decay constants and charge 
radii. In addition, it is also instructive and important 
to check the explicit dependence of these observables on 
the model’s parameter. 

Regarding these goals, we know from a previous 
study [17] that the value of the regulator mass, m^ = 0.6 
GeV, reproduces well all experimental data on the pion 
observables mentioned above. It is worthwhile to check 
whether = 0.6 GeV is also consistent with the kaon 
data, for which we compute the values of the decay con¬ 
stants and charge radii and compare them with their ex¬ 
perimental values. 

The calculated values of the observables listed in Ta¬ 
ble I show that mji = 0.6 GeV is also a suitable 
value for the kaon. Moreover, for mg = 0.44 GeV, 
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Decay Constant 

and 

Charge radius 


mu = 0.6 GeV; mg = 0.44 GeV 

Pion 


Kaon 

mu = md = 0.22 GeV 


mu = md — 0.25 GeV 


mu = md = 0.22 GeV 


mu = md = 0.25 GeV 

o o 

1 1 


93.12 MeV 

0.736 fm 


101.85 MeV 

0.670 fm 


110.81 MeV 

0.754 fm 


113.74 MeV 

0.687 fm 


TABLE I: Calculated decay constants and charge radius for two light-quark masses and corresponding experimental values. 
The experimental date with errors bar, are; = 92.42 ± 0.021 MeV, (r,r)™^ = 0.672 ± 0.008 fm, = 110.38 ± 0.1413 and 
= 0.560 ± 0.031 fm. Experimental date from [64, 70, 82]. 


mu = md = 0.22 GeV and mu = 0.6 GeV, the de¬ 
viation from the experimental values of f-K+{fK+) 
{r^+){{rK+)) are 0.76%(0.37%) and 8.7%(34.10%), while 
with = rud = 0.25 GeV, the mismatch is 9.3%(3.02%) 
and 0.35%(22.68%), respectively. 

For comparison, we recall that in case of a non- 
symmetric vertex function the ratio of decay constants, 
/xZ/tt, which measures 57/(3) flavor-symmetry break¬ 
ing, was found to be 1.363 [58]. Here, we read from the 
table that Jk/ f-w — 1.189, which is closer to the mea¬ 
sured value i.e., = 1.197 ± 0.002 ± 0.006 [82]. 

A graphic representation of the explicit dependence of 
the form factors, A®™, decay constants, /q- and charge 
radii (rp-) on the model’s parameters can be found in 
Figs. 1 to 8, in each of which one parameter is fixed while 
the other is varied. 

In Figs. 1 and 2, we plot the electromagnetic form fac¬ 
tors of the pion and kaon as a function of for various 
values of m/j and = rud and mg fixed. As can be 


seen, the elastic form factors in both figures are mono- 
tonically decreasing functions of with increasing hard¬ 
ness for larger values of mu. Fig. 1 also informs us that 
the available experimental pion data he in the interval 
of 0.1 GeV < < 1 GeV and from Fig. 2 we deduce 

that the kaon’s experimental form factor data are better 
reproduced for ^0.5 GeV which coincides with our 
privileged value, mu = 0.6 GeV. 

Similarly, in Figs. 3 and 4, the electromagnetic form 
factors of the pion and kaon are plotted as a function of q^ 
for different values of mu = md whereas mu = 0.6 GeV 
is fixed. One observes a likewise behavior of the elec¬ 
tromagnetic form factors, i.e. becomes harder 

for increasing values of mu = md- We note that for the 
symmetric vertex function with mji = 0.6 GeV, the light- 
quark mass should be in the range 0.15 < mq < 0.5 GeV, 
where the constituent quark mass, mq = 0.22 GeV, ap¬ 
pears to be the most favorable value to accommodate the 
experimental data. 

Next, the explicit dependence of the charge radii (r.^.) 



FIG. 1: The electromagnetic form factor of the pion as a 
function of space-like . The curves correspond to the 
different values of ttir with fixed quark masses: = 

md = 0.220 GeV. Experimental data: Ref. [66] (circle). 
Ref. [65] (square). Ref. [68] (diamonds). Ref. [69] (up trian¬ 
gle), Ref. [73] (down triangle). 



FIG. 2: The electromagnetic form factor of the kaon as a 
function of space-like q^. The curves correspond to the dif¬ 
ferent values of mu with fixed quark masses: mu = md = 
0.220 GeV, mg — 0.44 GeV. Experimental data: Ref. [83] 
(square). Ref. [64] (circle). 
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FIG. 3: The electromagnetic form factor of the pion as a 
function of space-like q^. The curves correspond to different 
values of nig = with regulator mass m_R = 0.60 GeV. Ex¬ 
perimental data: Ref. [66] (circle), Ref. [65] (square), Ref. [68] 
(diamonds). Ref. [69] (up triangle). Ref. [73] (down triangle). 



FIG. 5: Gharge radii (rg-) of the pion and kaon as a function 
of the constituent quark mass = md with nis = 0.44 GeV 
and fixed regulator mass mji — 0.6 GeV. The solid circle [70] 
and square [64] are the experimental values for the charge 
radii of the pion and kaon, respectively. 




FIG. 4: The electromagnetic form factor of the kaon as a 
function of space-like q^. The curves correspond to different 
values of niq = with fixed masses: rris = 0.44 GeV, mu = 
0.60 GeV. Experimental data: solid circle [64] (circle) and 
square [83]. 


FIG. 6: Gharge radii (rp-) of the pion and kaon as a function 
of the regulator mass mu with = md ~ 0.22 GeV and 
ms = 0.44 GeV fixed. The solid circle [70] and the square [64] 
are the experimental charge radii values of the pion and kaon, 
respectively. 


and (tk) on mu = md is depicted in Fig. 5 from which we 
deduce that the charge radius exhibits a nonlinear behav¬ 
ior and decreases strongly for large values of = rud, 
as expected. For instance, for m„ = md = 0.15 GeV, 
the size of the pion (kaon) charge radius is about 0.98 fm 
(^ 1 fm), whereas for = md = 0.5 GeV this size 
reduces to 0.43 fm (0.45 fm). A similar behavior can 


be seen in Fig. 6, where the charge radii are plotted 
against the regulator mass mji with mu = md = 0.22 
GeV and mg = 0.44 GeV. However, these effects are less 
pronounced for variations of mu than of = m^. 

Finally, the weak decay constants, /,r and Jk, are 
plotted in Fig. 7 and Fig. 8 as functions of = md 
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FIG. 7: The weak decay constants /q- of the pion (solid 
curve with rtiu = rrid = 0.22 GeV) and kaon (dashed curve: 
rria = 0.44 GeV) as a function of dotted curve: 

kaon decay constant as a function of rrig. In all curves mpt = 
0.6 GeV. The solid circle and square are the experimental 
decay constants of the pion and kaon, respectively [82]. 



FIG. 8: The weak decay constant, /q- , of the pion (solid 
curve) and kaon (dashed curve) as a function of the regulator 
mass mu with m^ = mj, = 0.22 GeV and = 0.44 GeV 
fixed. The circles and the square are the experimental values 
or the charge radii of the pion and kaon, respectively [82] . 


IV. CONCLUSIONS 


{niR = 0.6 GeV) and niR {niu = md = 0.22 GeV), where 
in both cases mj = 0.44 GeV. We observe that in contrast 
to the charge radii the decay constants are continuous in¬ 
creasing functions of rUg and niR. Moreover, the charge 
radii and the decay constants are more sensitive to the 
quark mass values, niq, than to mR. It is worthwhile to 
point out, as discussed in Ref. [58] for the non-syinmetric 
vertex, that the sensitivity of the kaon decay constant to 
the strange quark mass is very modest, whereas for the 
present symmetric vertex the mj dependence is quite sig¬ 
nificant. 

We stress that the regulator mass ttir = 0.6 yields the 
best fit to the experimental values of the decay constants 
and charge radii as discussed above. Furthermore, the 
decreasing (increasing) behavior of the charge radii (de¬ 
cay constants) with niq, as depicted in Figs. 5 and 7, 
satisfies Tarrach’s relation [84], i.e. (r-Q-.) ~ l/m, and 
/o- ^ 1 /(^ 0 -)- 


We reexamined the light-front approach to the light 
pseudoscalar mesons [16, 17, 58] by considering a sym¬ 
metric qq bound-state function [17]. In this framework 
and with the given symmetric ansatz for the bound-state 
vertex function, we calculated the charge radii, (r 7 r), and 
{’’’k), weak decay constants, /tt and and the electic 
form factors, F®™(g^) 7 r and F’^™{q^)K. 

To constrain our model parameters, namely, the quark 
masses m„, md and mg and the regulator mass mR, we 
first adjusted their values to repoduce the experimental 
weak decay constants. In doing so, we imposed the same 
regulator mass for both the pion and kaon and found 
that mR = 0.6 GeV is a suitable value to describe all 
experimental data on (’'ir(if)) ^^nd /^(ic) within 

the reasonable theoretical uncertainties. The numerical 
results also show this model significantly breaks down for 
TOij > 1 GeV, which was already demonstrated in the the 
case of a nonsymmetric vertex function [58]. 

Moreover, the explicit dependence of charge radii and 
decay constants on the quark masses, m„ = md and mj, 
not only satisfies Tarrach’s relation but also favors the 
range of mass values commonly chosen within the light- 
front model. In addition, by using the privileged val¬ 
ues of the model’s parameters m„ = md = 0.22 GeV, 
mg = 0.44 GeV and mR = 0.6 GeV, we find that the 
pion-to-kaon decay constant ratio is in excellent agree¬ 
ment with its experimental value, i.e. the mismatch is 
barely 0.67%. Lastly, the present numerical investigation 
suggests that these parameter values could also be used 
to study dynamical properties of other pseudoscalar and 
vector mesons or may apply to heavy-to-light transition 
form factors. 
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